Steady state Ab-initio Theory of Lasers with Injected Signals 
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We present an ab-initio treatment of steady-state lasing with injected signals that treats both 
multimode lasing and spatial hole burning, and describes the transition to injection locking or 
partial locking in the multimode case. The theory shows that spatial hole burning causes a shift in 
the frequency of free-running laser modes away from the injection frequency, in contrast to standard 
approaches. 
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Laser action in the presence of an injected signal is 
an extremely important topic both for research in non- 
linear dynamics and laser physics, and for applications 
of lasers. Under certain conditions, a laser mode can be 
locked to the injection frequency, allowing for stabiliza- 
tion and modulation of a "slave" laser based on control 
by a "master" laser. There has been great interest in such 
systems as a component in optical phase-locked loops 
for coherent detection [l| and for stabilization of novel 
laser systems at desired frequencies for optical commu- 
nications Q, El . However the conditions of stable locking 
are complicated by a fundamental principle of non-linear 
dynamics, that a driven, damped system of non-linear 
oscillators can exhibit complex dynamical behavior and 
even chaos when three or more independent, time- varying 
fields are coupled [4(. While the textbook theory of injec- 
tion locking, described by the Adler equation p, @|, does 
not describe such dynamical effects, most of the exten- 
sive research on injection locking in the past thirty years 
has focused on these issues [4] . However there is another 
aspect of injection locking which has not been treated at 
all in the literature, the space-dependent non-linear in- 
teractions of the injected and self-organized lasing fields. 
All previous treatments have immediately eliminated the 
spatial degrees of freedom through approximations which 
amount to treating the fields as uniform in space Q, 
hence completely neglecting the effect of spatial hole- 
burning, present in most real lasers, and making the the- 
oretical treatments qualitative at best. For many of the 
applications mentioned above, the systems are designed 
to avoid complex dynamics, and what is needed is a quan- 
titative theory of locking which takes into account the 
particular features of the laser cavities and gain media 
involved. 

To achieve such a theory, which includes the spatial 
variation of the injected signal and internal fields, it is 
critical to have a formulation which describes the open, 
non-hermitian nature of the system and the spatial varia- 
tion and interaction of the fields exactly. For multimode 
lasing without injection such a theory has been devel- 
oped in the past few years by Tureci, Stone and cowork- 
ers (il-[l2|. denoted by the acronym SALT (Steady-state 
Ab initio Laser Theory). The approach is valid for a wide 



range of laser systems, although it has been formulated 
mainly for those with atomic like gain media. It is a self- 
consistent theory in which only a few parameters of the 
gain medium are specified and all physical parameters of 
the lasing fields are calculated ab initio, e.g. the lasing 
thresholds, fields and frequencies. Here we show that this 
theory is generalizable to include arbitrary input fields 
and an arbitrary number of lasing modes, while remain- 
ing computationally tractable, so that injection locking 
can be quantitatively studied, even for cases with several 
lasing modes. This new theory, which we term Injection- 
SALT (I-SALT), assumes a free-running steady-state is 
achieved, and thus it is meant as an alternative to the 
Adler theory g, 0, [TM3 or other steady-state theories 
and not as a theory of the more complex dynam- 
ical states of lasers that can arise due to injected signals. 
We find below a novel effect of frequency repulsion from 
the injected signal which is not predicted by any pre- 
vious theory. I-SALT is shown to agree quantitatively 
with full time-dependent numerical solutions of the laser 
equations for a simple one-dimensional cavity with an 
injected signal, but because of its much greater compu- 
tational efficiency [19] . it will be feasible to apply it to 
realistic laser structures which are impractical to study 
using such direct simulation. 

SALT was formulated to determine directly the steady- 
state of the laser rate equations for an N-level atomic 
gain medium coupled to Maxwell's equations within an 
arbitrary cavity specified by its passive dielectric func- 
tion, e c (x), subject to a spatially- varying pump, D$(x) = 
doF(x), without performing time-integration to steady- 
state. It was shown that if a multi-mode steady-state 
exists, which requires stationary level populations, then 
all of the properties of the laser field are determined by 
a set of coupled time-independent wave equations with a 
saturating non-linearity that self-consistently couples the 
lasing modes. These equations are non-Hermitian at the 
linear level, and are subject to purely outgoing bound- 
ary conditions. They are solved efficiently by introducing 
a specific self-orthogonal basis set of threshold constant 
flux (TCF) states in which to expand the solutions, and 
then iteratively solving the resulting non-linear matrix 
equation for both fields and frequencies. We now sketch 



how this approach can be generalized to yield I- SALT. 

As in SALT, we assume the existence of a steady-state 
with stationary level populations; then the positive fre- 
quency components of the electric field inside the cavity 
for a given pump value, do, take the form 
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where the Nl lasing modes, have unknown spa- 

tial variation and unknown frequencies, c^, and there are 
Na amplified signals injected into the cavity, at given 
frequencies, cj a , and incoming amplitudes, B a , but with 
unknown amplitude and spatial variation, ^ a (x), within 
the cavity. All of the unknown quantities will be de- 
termined from the resulting I- SALT equations and their 
boundary conditions self-consistently. 

We have shown elsewhere that within the stationary 
population approximation (SPA) any N-level atomic laser 
system can be mapped onto an effective two-level model 
[19|, so here we consider an effective two- level Maxwell- 
Bloch model in which only the pump, do, and inversion 
relaxation rate appear, and the latter is absorbed into the 
units of the electric field [20| . Inserting the multi-periodic 
ansatz (pp) into the coupled equations for the field, polar- 
ization and atomic population inversion, and applying 
the SPA allows elimination of the latter two variables, 
leading to Nl + Na coupled non-linear wave equations, 
which can be written as three-dimensional vectorial equa- 
tions, but which here we will only consider in their scalar 
form, appropriate for the one-dimensional slab geometry 
studied below: 
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where Y v = 7 2 _/[(uv — uo a ) 2 + 7^] is the gain curve, 7^ 
is its linewidth, uo a is the atomic frequency of the lasing 
transition, k a = uj a /c is the wavevector, and a is an index 
running over both lasing states, /i, and amplified states, 
a. The equations for Sfr^x) are to be solved with purely 
outgoing boundary conditions, while those for ^ a (x) are 
to be solved with the boundary condition of fixed input 
amplitude at u a . 

We will solve these coupled equations by non-linear it- 
eration after expanding the solutions in a non-Hermitian 
basis set with the appropriate boundary conditions. For 
the lasing modes, this set is the same TCF states 
used in SALT, which satisfy, 

[V 2 + (e c (x) + rj n F(x)) k 2 ] u n (x; u) = (4) 
d x u n (x\u)\ x = L = iku n (L;uj), (5) 



where we refer to u n (x;uj) and rj n (oj) as the TCF eigen- 
vectors and eigenvalues respectively and for simplicity we 
have written the outgoing boundary condition for a one- 
sided cavity with a perfect mirror at the origin. Note that 
{r] n (uj)} is generically complex and can be thought of as 
the set of values of the gain medium susceptibility which 
lead to lasing, i.e. for which a solution for purely outgoing 
real wavevector exists. Thus the lasing thresholds in the 
absence of input signals are obtained by varying u in the 
TCF equation until a frequency, is found at which 
Vn(^fi) = d j±/(uJ^ - uJ a + *7Jl)- Tne {u n (x]UJ^)} then 
form an efficient basis set for finding the non-linear solu- 
tions above threshold because at the first lasing thresh- 
old, the lasing mode is just a single TCF state, and above 
threshold only a small number of TCFs are needed to con- 
verge to the non-linear solution of the SALT equations. 

The amplified modes, however must be treated 
differently from the lasing modes since they have a fixed 
incoming signal amplitude and fixed frequency. To repre- 
sent these modes we require terms with an incoming com- 
ponent in addition to the outgoing TCF expansion terms, 
which we do conveniently by solving the same TCF equa- 
tion inside the cavity with a purely incoming boundary 
condition, 

= ^2 a n a) ^nO; Ua) + ^ b m )v m(x] ^a), (6) 
n m 

where the states v m (x;uj) and associated eigenvalues /3 m 
are given by 

[V 2 + (e c (x) + P m F{x)) k 2 ] v m (x; u>) = (7) 
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and thus represent states that are purely incoming. The 
incoming and outgoing TCF states are not power orthog- 
onal, but they do satisfy a self-orthogonality condition 
between themselves [2l| . Either the incoming or outgoing 
TCF states represent a complete basis for fields within 
the cavity at uj ai but the incoming terms are needed to 
represent the input boundary condition. Because they 
are purely incoming, they do not contribute directly to 
the emitted fields, but they correctly represent the full 
spatial hole-burning and gain competition effects of the 
amplified input. 

For amplified modes, we can easily write the incoming 
boundary condition for a one-sided slab cavity of length 
L as 



-ik a L 
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(9) 



where B a is the given incoming field amplitude at fre- 
quency uj a . This single equation vastly underdetermines 
the coefficients b^m in the sum, so that the choice is 
based on convenience. This freedom arises from the over- 
completeness of using both {u n } and {v m } to represent 



the internal fields. Hence the coefficients depend 
strongly on the choice of the bffl • A natural choice is 
to take only a single term, vo(x;u a ), which corresponds 
to the dominant component of the outgoing TCF state 
for the nearest lasing mode. This is allowed for a cavity 
with a single input channel, as in the one-sided slab ge- 
ometry we are considering here; in general one needs a 
minimum of M independent incoming states to represent 
an arbitrary input for an M-channel cavity, and these can 
be chosen again to be similar in character to the nearest 
lasing mode in order to optimize the calculation. 

Once a representation of the input field is chosen, one 
can insert the expansions ((6]) for the amplified modes and 
the analogous expansion involving only outgoing TCFs 
for the lasing modes into the fundamental Eqs. ((2]) and 
(j3j) and use the self-orthogonality relations of the u n , v m 
to find coupled non-linear matrix equations for the coeffi- 
cient a^n\ bm ^ which determine their solutions. For 
the lasing modes, one finds 
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where 7 M 7±/(^ 



*7_l)- This is identical to the 



lasing equations of SALT except for the presence of the 
amplified mode intensities in the non-linear hole-burning 
denominator. 

In a similar manner the coupled equations for the am- 
plified modes can be determined, and they take the form 

m-. 

= E T ln )a i a) + E + V^) (12) 

(13) 



V£ m J = — - —ui(L] uj a )v m (L; uj a ). 



The operator in ([T2]) takes precisely the same form 

as except that u n {x) — > v m (x). 

Together Eqs. (10-13) define I-SALT. In the regime in 
which the SPA holds, they provide essentially exact solu- 
tions of the full coupled wave equations for amplification 
and injection locking, which has not been achieved pre- 
viously, to our knowledge. The method is ab initio, as 
in SALT, with no prior assumptions about the number, 
spatial form or frequencies of the lasing modes. The phys- 
ical picture of injection locking in this regime is one of 
cross-saturation, and not synchronization of non-linear 
oscillators, as it is in the Adler theory. Lasing modes 
correspond to poles of the non-linear scattering matrix 
on the real axis; amplified inputs are simply additional 
scattered waves which also deplete the gain. If the in- 
put signal becomes too strong, and is sufficiently near in 
frequency to the lasing mode, then the lasing mode has 
insufficient gain and falls below threshold (the pole falls 



below the real axis [22]), leaving only the amplified signal 
output. The output is "locked" to the input frequency, 
but not by pulling the lasing mode over to u; a , but rather 
by turning it off. This scenario is appropriate when all 
the relevant frequencies or rates in the problem are large 
compared to the relaxation rate of the gain, so that co- 
herent interactions of the fields are averaged away. This 
is a common regime for solid state lasers with atomic 
like gain media and essentially the same conditions hold 
in the stable locking regime of semiconductor lasers. 

To test the results of I-SALT, in Fig. [T] we compare 
them to the exact numerical solutions based on FDTD for 
a simple slab cavity with an injected signal (schematic) 
under the standard conditions of a single lasing mode in 
the absence of the signal. Locking of the output signal 
is shown in the top panel based on I-SALT, in excellent 
agreement with FDTD data, with no adjustable param- 
eters. The value of the output signal leading to locking 
found by I-SALT is approximately the same, but larger 
than predicted by the Adler theory. In Adler theory one 
finds that the energy output from the gain medium at the 
locking threshold and in the free-running state with no 
input are exactly equal, and typically this is interpreted 
as due to "energy conservation" . But this is not quanti- 
tatively correct in a more complete theory, because, due 
to spatial hole-burning, it is possible for two output fields 
to extract more energy from the gain medium in steady- 
state than a single one, and indeed we find that the am- 
plified output at locking in Fig. 1 is ~ 18% larger than 
the original free-running output. 

Moreover, the frequency behavior of the lasing mode 
prior to locking is qualitatively opposite to the Adler the- 
ory; the frequency is repelled from the input frequency 
instead of being attracted to it. The repulsion here is dis- 
tinct from that observed by Murakami et al. [15| . arising 
from the Henry a-factor, which is negligible in atomic 
gain media. This frequency repulsion can be explained 
by the combined effects of mode competition and spa- 
tial hole-burning. As the incident signal is imposed and 
depletes the gain, the standing wave of the laser field 
shifts away from the frequency /wavelength of the inci- 
dent standing wave in order to better extract energy from 
the regions of the cavity where the gain is not being sat- 
urated by the incident signal. To confirm this interpre- 
tation we replace the space dependent gain saturation 
denominator with its constant spatial average in the I- 
SALT equations and find no substantial frequency shift. 

As noted, the Adler theory describes locking driven 
by phase synchronization of the input and free-running 
fields. Since the threshold input intensity for locking de- 
creases to zero as the input frequency approaches the 
free-running frequency, the threshold condition can be 
expressed as a "locking range", the frequency range Acj, 
over which the laser is locked for a given input inten- 
sity. In the Adler theory one finds Acj = ^ c ^/\Ba\^/I^^ 
where j c is the cavity decay rate, \B a \ 2 is the intensity of 
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FIG. 1. (Top panel) Injection locking transition as described 
in the text. Solid lines are output intensities calculated from 
I-SALT; blue is lasing output, red is amplified output at sig- 
nal frequency, dot-dashed black is total output. Triangles are 
the same quantities from FDTD for the same dielectric slab 
laser with uo a = 40, 7_l = 4 (schematic inset). To the left of 
the black dashed vertical line the laser is free-running with 
no input signal and output at 0Jij ree = 40.714; to the right 
an input signal at uji n — 40.4 is applied and increased (see 
top inset), eventually driving the lasing signal to zero at the 
locking transition (orange dashed line). (Bottom panel) Fre- 
quency variation of the first lasing mode: blue line is from 
I-SALT and blue triangles from FDTD. The green line shows 
the prediction of the Adler equations, the green circle show- 
ing where locking is predicted. Again, the orange dashed line 
shows locking threshold from I-SALT, frequencies beyond this 
point are taken as the real part of the location of the pole of 
the scattering matrix. As the locking transition is approached 
the lasing frequency moves away from its free-running value, 
due to spatial hole-burning in contrast to the behavior ex- 
pected from the Adler equation [5]. Blue dashed lines shows 
I-SALT calculation with spatially averaged hole-burning. The 
inset shows a plot of the phase shift between input and output 
signals of an injection locked dielectric slab cavity at a fixed 
input intensity. I-SALT (red) and FDTD simulations (red tri- 
angles) are seen to have close agreement and a smaller locking 
range than the Adler prediction (green). Values quoted are 
given in units of c/L. 



the input signal, and i^o is the intensity of the free run- 
ning lasing signal in the absence of the input [5| . Within 
this locking range there is a fixed phase relationship be- 
tween the input signal and the locked output which varies 
as Acj) = arcsin^^o — co^/Acj]. 

The same quantity can be calculated in I-SALT and is 
compared to the Adler prediction for the same slab cav- 
ity in the inset of Fig. [TJ The phase shift variation found 
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FIG. 2. Partial locking transition as described in the text 
using a similar pumping and input ramping scheme as Fig. 
[U Solid lines are output intensities calculated from I-SALT; 
blue and cyan are lasing output, red is amplified output at 
signal frequency, uj = 20.3, dot-dashed black is total output. 
Triangles are the same quantities from FDTD for the same 
dielectric slab laser with uo a — 20.5, 7^ = 3 (schematic inset). 
First the nearby lasing mode locks to the injected signal (or- 
ange line), then the distant lasing mode locks (purple line). 
Values quoted are given in units of c/L. 

from I-SALT is completely different from the Adler pre- 
diction and in good agreement with FDTD. Note that 
I-SALT also finds a slightly smaller locking range than 
Adler theory, a reflection of the fact that the input in- 
tensity for locking is larger than that predicted by Adler 
theory. 

Beyond yielding a better treatment of the locking tran- 
sition of a single mode laser, as seen from the generality 
of Eqs. (10-13), I-SALT is able to treat simultaneously 
multimode lasing with multiple inputs. An example of 
this with two mode lasing and a single input signal is 
shown in Fig. [2] Here a similar slab laser is pumped 
above its second lasing threshold and a signal is injected 
closer to the frequency of the first lasing mode. Because 
of its stronger interaction with the first mode (blue), the 
signal is able to lock that mode, while the second mode 
(cyan) remains active at a similar frequency to its free- 
running value, though shifted away from the injected fre- 
quency in the same manner as described before. Thus 
we have a "partially-locked" state of a type not treated 
in any previous theories. Again, excellent agreement is 
found between I-SALT and FDTD, but with substan- 
tially less computational effort, I-SALT takes an hour 
of computational time on modern CPUs to generate the 
curves seen in Fig. [2] while FDTD requires 168 days. 
In general, exact FDTD studies of multimode lasing are 
computationally very demanding and could not be per- 
formed in more realistic structures, whereas SALT has 
been shown to be computationally tractable in complex 
two-dimensional structures such as photonic crystals [23| 
and random lasers [To[, and 3D vectorial codes are cur- 
rently under development and look promising. Since I- 
SALT is essentially of the same degree of computational 
complexity as SALT, it can be used as a tool to study 



frequency control of both single and multimode lasing, 
and opens up new possibilities for investigating injection 
locking systems. 
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